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Abstract: The D-instanton corrected hypermultiplet moduli space of type II string theory com- 
pactified on a Calabi-Yau threefold is known in the type IIA picture to be determined in terms of 
the generalized Donaldson-Thomas invariants, through a twistorial construction. At the same time, 
in the mirror type IIB picture, and in the limit where only D3-Dl-D(-l)-instanton corrections are 
retained, it should carry an isometric action of the S-duality group SL(2, Z). We prove that this is 
the case in the one-instanton approximation, by constructing a holomorphic action of SL(2, Z) on 
the linearized twistor space. Using the modular invariance of the D4-D2-D0 black hole partition 
function, we show that the standard Darboux coordinates in twistor space have modular anomalies 
controlled by period integrals of a Siegel-Narain theta series, which can be canceled by a contact 
transformation generated by a holomorphic mock theta series. 
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1. Introduction 



In type II string vacua with N = 2 supersymmetry in four dimensions, the moduli space spanned 
by massless scalar fields famously factorizes as the product M. v x M. u of the vector multiplet 
(VM) and hypermultiplet (HM) moduli spaces, respectively 0,0. A complete understanding of 
the former was achieved by the end of the last millenium, and has led to key developments in black 
hole physics and algebraic geometry. By contrast, until recently our understanding of the latter 
was confined to the tree-level |3], |, |5|] and one-loop || 0, ^ |^, [Kj approximations to the metric, 
together with some general expectations about the form of instanton corrections [II], |T2]]. While 
it has been clear to many that an exact solution for the quaternion-Kahler (QK) metric on Ain 
would have considerable impact both for physics and geometry, progress was hindered mainly by 
the lack of a convenient parametrization of this class of metrics. 

The situation has improved dramatically in recent years, as projective superspace and twistorial 
techniques [O, [L4l [L5l [TBI O, llSL [19] were brought to bear on this problem. The key idea behind 



these methods is that the QK metric on A4h is encoded in the complex contact structure on the 
twistor space Z, a P 1 bundle over M. H [JTE| . The latter is in turn specified by a set of holomorphic 
generating functions for complex contact transformations between local Darboux coordinate sys- 
tems, subject to certain global consistency and reality requirements [F2DJ. These gluing conditions 
can be converted into a system of integral equations, which can often be solved by iteration. The 
QK metric on jW^ is then obtained by expanding the resulting Darboux coordinates around any 
fixed section of Z, similar to the usual twistorial construction of HK manifolds |15|, |21 |. 



Using these techniques, combining earlier results on the twistorial description of the pertur- 
bative metric and of Dl-D(-l)-instantons [22, 23, ^ |25| and taking inspiration from a similar 
construction in the context of gauge theories [23], a general construction of the D-instanton cor- 
rected HM moduli space was laid out in [|27|, [28|], in terms of the generalized Donaldson-Thomas 
(DT) invariants f2(7; z) which count D-instantons. For reasons that will become clear below, we 
refer to the construction of |E7|, as the 'type IIA' construction. In type IIA string theory com- 
pactified on a Calabi-Yau (CY) threefold 2), D-instantons come from D2-branes wrapping special 
Lagrangian cycles, while in type IIB string theory compactified on the mirror threefold 2J, they 
come from D5-D3-Dl-D(-l)-branes wrapping even-dimensional cycles, or more generally coherent 
sheaves on 2J. Although the DT invariants typically depend on the complex structure moduli (in 
type IIA) or Kahler moduli (in type IIB), they are piecewise constant away from walls of marginal 
stability, and their discontinuity is such that the resulting contact structure on Z (and hence, the 
QK metric on M.h) is smooth PB[ . In fact, the type IIA construction of [^7], is essentially 
dictated by consistency with wall-crossing, and is isomorphic to the gauge theory construction of 
2"6| ] under the QK/HK correspondence [2_5|, |3D . 



In addition to these D-instantons, one also expects contributions from NS5-brane instantons 
wrapped on the CY threefold. Those effects are in principle determined from the known D5- 
instanton corrections, by requiring that the QK metric should carry an isometric action of SL(2, Z), 
originating from S-duality in ten-dimensional type IIB string theory. NS5-brane corrections were 
constructed at linear order in [3JJ (see also JS2J, |33J, |3jj[), but an understanding of these effects at 
non-linear level is still missing. 

In this work, we focus on a subset of the allowed instanton corrections, namely those corre- 
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sponding to D3-Dl-D(-l)-instantons on the type IIB side. Although these effects are exponentially 
suppressed compared to Dl-D(-l)-instantons in the large volume limit, they are still exponentially 
larger than D5 and NS5-brane instantons. Thus, a natural question is whether the twistorial con- 



struction of p7| , Pq] produces a modular invariant quaternion-Kahler metric, in the large volume 
limit where D5-NS5- instantons are ignored but the effects of D3-Dl-D(-l)-instantons with non-zero 
D3-brane charge are still retained. 1 In this work, we shall answer this question in the affirmative, 
at least in the one-instanton approximation. As will become clear shortly, even in this simplified 
setting, modular invariance is achieved in a very non-trivial way, and depends on the sophisti- 
cated machinery of Eichler integrals and indefinite theta series. We expect an even richer structure 
beyond the one-instanton approximation, but we shall hardly touch upon it in this work. 

The reason why this problem is non-trivial is that modular invariance is far from manifest 



in the 'type IIA' twistorial construction |27, Indeed, when lifted to a holomorphic action on 
twistor space, the action of S-duality involves a SU (2) rotation along the fiber, and therefore relates 
Darboux coordinate systems on different patches. The only exception are the patches around the 
points t = ±i, where t is the usual stereographic coordinate on the fiber, which are mapped onto 
themselves under S-duality. But even in these patches the Darboux coordinates need only be 
modular covariant up to a local complex contact transformation. Finally, a technical but serious 
difficulty is that the natural type IIA coordinates parametrizing the twistor fiber need not have 
simple transformation properties under S-duality, and it is a priori unclear how to construct a basis 
of type IIB fields which would transform covariantly. Put differently, there is no general method to 
construct the quantum mirror map (see however for recent progress in a closely related set-up) . 

Fortunately, we find that this picture drastically simplifies in the large volume limit V — > oo, 
provided one simultaneously scales z ~ (t ± i) — > keeping the product zV 1 ^ 3 fixed. In this 
case all BPS rays (the contours on P 1 associated to D-instantons [^6], |27fl ) coincide with the real 
axis in the complex z plane, shifted slightly above or below the origin z = depending on the 
Dl-brane charge. The corrections to Darboux coordinates are then given by integrals of certain 
Siegel-Narain theta series along these contours, corresponding to the sum over Dl-brane charges 
at fixed D3-brane charge. These Gaussian theta series are far simpler than the non-Gaussian sums 
which would otherwise arise at finite volume. 

Another simplification of the large volume limit is that the moduli-dependent DT invariants can 
be traded for another set of invariants, which we call the MSW invariants, which are independent of 
the Kahler moduli in type IIB (or complex structure moduli in type IIA), and have simple modular 
properties. The idea is that D3-Dl-D(-l)-instantons are governed by the same DT invariants which 
count D4-D2-D0 black holes in type IIA string theory compactified on the same CY threefold 2). At 
large volume, the latter can be described by micro-states of the MSW superconformal field theory 
(SCFT) ||37| , and bound states thereof |39| . We define the multi-instanton expansion in the 
large volume limit in such a way that only elementary MSW states contribute in the one-instanton 
approximation. With this definition, the Gaussian theta series alluded to in the previous paragraph 
reconstruct the elliptic genus of the SCFT (or modular derivatives thereof), and the corrections 
to Darboux coordinates are now identified as Eichler integral of the MSW elliptic genus. For 



1 For zero D3-brane charges, i.e. in the limit where only Dl-D(-l)-instantons are retained, the QK metric and 
twistor space are known to be modular invariant [36fl , although this is not apparent in the usual type IIA 

formulation. 
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illustration, we display the result for the correction to the Darboux coordinates £ conjugate to the 
Dl-brane charges q, evaluated for simplicity at the point z = on the twistor fiber, 



4?r ..rr^/A a/i(w-t) 



/xeA*/A 



where S c i is the classical D3-instanton action ( [4.11 ). It involves the Siegel-Narain type theta series 
d p ,fj,(w, f , t, b, c) ( |A.16| ) of weight (|, ^f^), where 6 2 = ^(2)), analytically continued away from the 
slice w — t, and the partition function of the MSW invariants h p /1 (T), which is a vector- valued 
modular form of weie ht (-|-1,0). As a consequence, the Darboux coordinates around t = ±i are 
not modular covariant, but transform with a specific modular anomaly given by the period integral 
of the MSW elliptic genus. In the example above, this implies 



+ «(0) + f e~™* £ h pA r) r . (1.2) 



/xeA*/A 



The final step in the proof is to show that the modular anomalies in the Darboux coordinates 
can be absorbed by a complex contact transformation, which will then ensure that the complex 
contact structure on Z transforms covariantly under S-duality. For this purpose, we observe that 
the Eichler integral which determines the corrections to the Darboux coordinates is in fact the 



modular completion of the holomorphic mock theta series introduced in [40 



9 PiM (r, t, t', b, c) = (-V k ' P ^ k + 6 ) • *) " s ^ k + 6 ) ' O] e- m ^ 2+2 ^ k+ 2 b \ 

(1.3) 

where t are physical Kahler moduli and t' is an arbitrary reference point on the boundary of the 
Kahler cone. In an amusing role reversal, we find that this holomorphic mock theta series provides 
the generating function of a complex contact transformation which cancels the modular anomaly 
in the Darboux coordinates. This establishes the modular invariance of the QK metric in the large 
volume, one-instanton approximation. 

Clearly, an important challenge for the future is to establish modularity at the non-linear level. 
One possible strategy would be to recast the type IIA twistorial construction into the manifestly 
SL(2, Z)-invariant type IIB formalism presented in ||36|| . As a first step in this direction, one may 
try to use the Poincare series representation (also known as Farey tail 33 1) of the generating 



functions of the MSW invariants to represent the section of H l (Z, 0(2)) describing D3-instanton 
corrections at linear order as a sum H = J2 m n Gm,™ suc h that S-duality acts by permuting these 
contributions. Unfortunately, this manipulation is formal since Poincare series with negative weight 
are divergent without a suitable regularization. Even though this naive idea fails, it nevertheless 
provides very useful guidance for determining the quantum mirror map, as we explain in Appendices 
B and C. 

The remainder of this work is organized as follows. In §[| we review the general properties of 
the HM moduli space in type IIB string theory compactified on a CY threefold, with particular 
emphasis on the discrete symmetries under S-duality and large gauge transformations, and recall 
the type IIA twistorial construction of D-instantons. In §13], we recall the geometric description of 
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D3-brane instantons in terms of coherent sheaves, the definition of the generalized DT invariants 
which count them and their relation to the MSW invariants, which are independent of the moduli 
and have simple modular properties. We define a multi-instanton expansion which is best suited 
for studying modularity, and give a heuristic explanation of modular invariance in twistor space. 
In we analyze the instanton corrections to the Darboux coordinates in the large volume, one- 
instanton limit, zooming in near the S-duality invariant points t = ±i. We express these instanton 
corrections in terms of Eichler integrals of Siegel-Narain theta series, and construct their modular 
invariant completion. Some open problems are discussed in §|[ The relevant properties of indefinite 
theta series and Eichler integrals are reviewed in Appendix [A], while Appendix [FJ contains details 
about the derivation of instanton corrections and quantum mirror map in the large volume limit. 
Finally, Appendix |CJ gives a preliminary attempt to recast the type IIA twistorial construction of 
D3- instantons in the type IIB framework of |3TSfl . 



2. HM moduli space in type IIB CY vacua 

In this section, we briefly recall some basic facts about the hypermultiplet moduli space in type IIB 



CY vacua. We mainly follow p7|| , incorporating further improvements introduced in [ 28J, |3l], p0|| . 



as well as some new insights on the action of S-duality on the tree level geometry. 
2.1 Tree-level metric 

Recall that the hypermultiplet moduli space in type IIB string theory compactified on a CY three- 
fold 2) is a quaternion-Kahler manifold M.h of dimension 4&2+4, where b 2 = &2(?J), which describes 
the dynamics of 

1. the ten-dimensional dilaton g s ; 

2. the Kahler moduli b a + it a = f J (a = 1, . . . , 62) where J = B + i J is the complexified 
Kahler form on 2J and 7 a is a basis of H 2 (2J , Z) ; 

3. the Ramond-Ramond (RR) scalars c°,c a ,c a ,co, corresponding to periods of the RR 0-form, 
2-form, 4- form and 6-form on a basis of H even (%), Z); 

4. the NS axion if), dual to the NS 2-form B in four dimensions. 

The ten-dimensional string coupling r 2 = l/g s and the RR axion T\ = c° combine into the ten- 
dimensional axio-dilaton field r = T\ + ir 2 . The resulting field basis r, t a , b a , c a , c a , c , if) is adapted 
to the action of S-duality in ten dimensions, as will become apparent below. 

At tree level in type IIB string perturbation theory, the metric on is obtained from the 
moduli space Ais/c of complexified Kahler deformations (which also describes the vector multiplet 
moduli space in Type IIA string theory compactified on the same CY threefold 2J) via the c-map 
construction |]3], FJj. In the large volume limit, the special Kahler manifold Aisic is characterized by 
the holomorphic prepotential 

F cl (X A ) = -n abc + \A A ^X^ (2.1) 
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where X A are homogeneous complex coordinates on Msk such that X A /X° = z A (with z° = 1), 
K a b c is the triple intersection product in #4(2), Z), and A\y> is a constant, real symmetric matrix 
which satisfies the quantization conditions PT|] 2 

A 00 eZ, A Qa = ^ + Z, A ab p»-^ abc p b p c eZ forVp a GZ. (2.2) 

The c-map construction || produces the QK metric 
dr 2 



,7,2 

aS M„ - r 2 



+ Ads 2 SK - 1 (dC A - A/" AS dC S ) Im A/" AA ' (d( A ' - AWd^' 



+ 16^( da + 4dCA - CAd4 ) ' 
where ds^ is the metric on AisK with Kahler potential K = — \og[i(X A F A — X A F^)], 



(2.3) 



AAA' = ^AA' + 2i —=■, (2.4) 

X Imrss'-^ 

is the Weil period matrix, r A £ = d x ^d x ^F cX (X) the Griffiths period matrix, and 

r = %V, V = \\ J NJ NJ = \ K abc t a t b t c . (2.5) 
2 6 6 

Here V is the volume of the threefold 2) in string units, and the variable r is related to the four- 
dimensional string coupling by r = l/g\. The variables z a ,( A ,(A,& appearing in ( |2.3|) are the 
natural variables in type IIA string theory compactified on the threefold 2) mirror to 2J, and are 



related to the type IIB fields r 1; t a , b a , c a , c a , c , ip by the 'classical mirror map' [44 



z a = b a + it a , C° = n , C = ~(c a - nb a ) , 
1 1 

Ca = Ca + ~ ^abc b\c C - 7i& C ) , Co = C ~ g ^abc '(c° - T X b C ) , ( 2 .6) 

o = -2(4> + ~ Tl co) + ~c a {c a - nb a ) - \ K ahc b a c\c c - n b c ) . 

Here the primed coordinates Ca = Ca — ^4asC A are the RR fields in a non-integer homology basis 
on the type IIA side, which is best suited for comparison with type IIB f3~Ifl . The electric charges 
q' A in this basis are related to the electric charges g A in the integer basis by q' A = q\ — v4 As p s , and 
therefore satisfy the quantization properties 

p° 1 1 

q'a e Z - ^ C 2,a - ^abcV h V\ % G Z - ^jP^a- (2.7) 

The prepotential in the non-integer basis is given by ( |2.1| ) with Aas — 0, which we henceforth 
denote by F ,cl (X A ). 



2 For simplicity we set Aoo = in this paper, as this can always be achieved by an integer symplectic transfor- 
mation. The constants C2, Q are the components of the second Chern class of 2J along the basis of #4(2), Z). 
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2.2 Continuous symmetries at tree-level 

The tree-level metric ( |2.3|) is invariant under a large group of continuous isometries which is gen- 
erated by three subgroups: 

• continuous shifts of the RR and NS axions 

: (C A , Ca, <t) H- (C A + V A , Ca + VA, <r + 2k - r/ A C A + tj a Ca) , (2.8) 

• monodromies around the large volume point 



Mea : b a ^b a + e a , C a ^C° + e°C°, C ^ C " KafcC^ - ^ abc eVC° 
Co ^ Co - C> a + \ + ~ K abc e a e b e%° ; 

continuous SX(2,R) transformations given in the type IIB field basis by 



(2.9) 





+ b 


T !-)■ 


+ d' 


cr 










b a 





t a i — ^ t a |cr + d| , c a i->- c, 



a • 



V J V-6 a J U 



(2.10) 



with ad-bc=l M, 44 



The algebra generated by these isometries is the semi-direct product SL(2, R) x AT, where 

AT = j\T(i) © AT(2) ® jV(3) 

is a nilpotent algebra of dimension 3&2 + 2, satisfying 

[iV (1) ,iV (1) ] c iV (2) , [N {p) ,N {q) ] = ifp + g>3, (2.11) 

with the generators in , , transforming as b 2 doublets, b 2 singlets and one doublet under 
SL(2,M), respectively. The relation of this algebra decomposition to the symmetries presented 
above is the following: 

the group elements Tj^}^ obtained by exponentiating (with a suitable admixture of 



higher level generators) consist of the monodromies ( |2.9| ) and the Heisenberg shifts (|2.8|) with 
non- vanishing r) a ; 

• the group elements obtained by exponentiating correspond to the Heisenberg shifts 
with non- vanishing fj a ; 

• the group elements ^ associated to coincide with Heisenberg shifts with non- vanishing 
fj and k; 

• unlike the other Heisenberg shifts which belong to N, the Heisenberg shift (° h-> (° + 77° 
belongs to SX(2,R) and in fact coincides with the SX(2,R) transformation r 1— > t + rf . 

In the presence of quantum corrections, all these continuous isometries are lifted, but a discrete 
subgroup of these symmetries is conjectured to remain unbroken, as will be discussed in §EO. 
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2.3 Twistorial description of the tree-level metric 

In studying instanton corrections to hypermultiplet moduli spaces, it has proven very useful to 
encode the metric on a QK manifold M. in terms of the complex contact structure on its twistor 
space Z, the total space of the P 1 -bundle over M. twisted with the projectivized SU (2) connection 
[46(1 for a review of this approach). This contact structure is represented by a (twisted) 



sec 



holomorphic one-form X, which locally can always be expressed in terms of complex Darboux 
coordinates as 

#M = d«W + il ] d$ } , (2.12) 

where the index [i] labels the patches Ui of an open covering 3 of P 1 . The contact structure on Z 
(hence, the QK metric on Ai) is then encoded in holomorphic generating functions H^(£&,(ffl, a^) 
for the contact transformations between Darboux coordinate systems on the overlap Ui (lUj. For 
QK manifolds with one quaternionic isometry, which in our case is realized as constant shifts of 
the NS axion a, the transition functions should be independent of a^. In this case, the 
corresponding contact transformation takes the following form 

LJJ 1 J ?A S M /o 1 o\ 

and reduces to a complex symplectomorphism in the subspace spanned by (£ A , £a) variables. In 
addition to the set of transition functions one must also specify a real constant c a , known as 
the anomalous dimension, which characterizes the singular behavior of a near the poles of P 1 . The 
QK metric is then obtained by 'parametrizing the twistor lines', i.e. solving the gluing conditions 
( p,13| ) for the Darboux coordinates as functions of coordinates on the base Ai and the P 1 coordinate 
t, and expanding the contact one- form in the vicinity of any fixed point on P 1 . A useful construct 
in this procedure is the contact potential e*, which determines a Kahler potential on Z (see [^0] 
for more details on this procedure). 

As explained in [^, |27|, the QK metric Q2.3| ) can be cast in this twistorial framework by 
choosing a covering of P 1 consisting of two patches U + , U- around the north and south poles, t — 
and t = oo, and a third patch Uq which covers the equator. The transition functions between 
complex Darboux coordinates on each patch are given by 

H [+o] = f c \£ a ), H [ -° ] = F cl (£ A ), (2.14) 
whereas the corresponding Darboux coordinates in the patch Uq read as B7I, PD 



a 

-cl 



fA = Cv + i (t-'Ft-tF*), (2.15) 
o + f {t- x W cl - t W cX ) , 



where W denotes the 'superpotential' 

W c \z)=F£{z)C K -z K U (2.16) 



3 In principle, since Z is a non-trivial P^bundle over A4, one should be using an open covering of Z. However, 
for QK manifolds with one quaternionic isometry, as is the case in this paper, the QK/HK correspondence gives a 
natural trivialization of an open subset of Z, which allows to view Z as a .M-bundle over P 1 . 
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and we have traded the Darboux coordinate a for a = — 2a — £a£ a . Moreover, the contact potential 
e* turns out to coincide with the four- dimensional string coupling r given in ( |2.5| ). 

A very useful property of this twistorial framework is that isometries of M. can always be lifted 
to a holomorphic contact transformation on the twistor space Z, by a suitable choice of action on 
the P 1 fiber. Thus, the nilpotent generators ( |2.8| ) and ( |2.9| ) (supplemented with the trivial action 
on P 1 ) lift to the holomorphic contact transformations 

V.fcv,*) : (£ A >£~A,a) (e A + ^ A , U + VA, a + 2K-^ A + V A U), (2.17) 



i'o Co - iae a + I K a U a e b e c + ~ K abc e a e b e c e , a^a, 
2 o 



(2.18) 



respectively, where we defined the Darboux coordinates in the non-integer symplectic basis in the 
same way as in (p.6[), 

? A = a - A AJ ,e, a' = a- 1 - A^ee- (2.19) 



Similarly, the isometric SL(2,M.) action ( 2.10 ) lifts to the following holomorphic action on Z 
(here expressed in terms of the Darboux coordinates in the patch Uq): 

l'o\ ^ (d -c\ /g\ 1 6 / cV(c^ + d) ^ (2 -- U) 

«7 \-b a J \a'J + 6 abcK * 5 V-[c 2 «° + 6) + 2c]/(c£° + d) 2 

This action agrees with Q2.10| ) provided the latter is supplemented with a SU(2) rotation along the 
fiber, 

cf + d , . 

m, ^z, 2.21 

where the coordinate z is related to the coordinate t appearing in ( |2.15| ) by the Cayley transfor- 
mation: 

t + i + - 1 + z (ooo\ 

z = r , t = -l . (2.22) 

t — i 1 — z 

Under this action, the complex contact one-form (|2.12| ) transforms by an overall holomorphic factor 
i — ^ Af[°y(c^ + cf), leaving the complex contact structure invariant, while the contact potential 
e* transforms with modular weight (—§,—§), 

e*^^— . (2.23) 
|cr + d\ 

It will be key to note that the points z = 0, oo (corresponding to t = ±i) are invariant under 
SL(2, R.). As a result, the leading Taylor coefficients of the Darboux coordinates around these 
points have simple modular transformations, e.g. £ a = rb a — c a + O(z) transforms with modular 
weight (—1, 0). 
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2.4 Quantum corrections and S-duality 

Away from the large volume, weak coupling limit, the continuous SL(2, R) isometric action is broken 
by several types of quantum corrections. Firstly, due to worldsheet instantons the holomorphic 
prepotential ( |2.1| ) entering the c-map metric ( |2.3| ) is corrected into F = F cl + _p w s - where 

C(3)(X ) 2 (X ) 2 ^ 
F {X ] = X * 2(2*1)3 ~ (2^ ^ ^ Ll H C J' (2 ' 24) 

The first term, proportional to the Euler number x^ of 2) , corresponds to a perturbative correction 
in the a'-expansion around large volume. The second term corresponds to a sum of worldsheet 
instantons, labeled by their effective 4 homology class k a ^ a e H^ii), Z), and weighted by the genus 
zero Gopakumar-Vafa invariants nj^ G Z. These instantons contribute through the tri-logarithm 
function Li^(x) = Ylm=i m ~ 3xm '■> w hich takes into account multi-covering effects. Note that the 
two terms in ( |2.24| ) may be combined by including the zero class k a = in the sum and setting 



71, 



(0) 




— Xm/2. In the presence of these worldsheet instantons, the twistor space is still described 
by three patches related by the transition functions ( J2.14TJ , where F cl (£ A ) is replaced by _F(£ A ). 

Secondly, the HM moduli space receives corrections controlled by the string coupling g s . At 
the perturbative level, it is believed 0, §, ^Of that the only corrections occurs at one-loop, and is 



proportional to the Euler number x^ H 0- Its effect on the metric was explicitly calculated in 
H |l0| . At the level of the twistor space, the one- loop correction is simply incorporated by allowing 
a logarithmic singularity in the Darboux coordinate a with prescribed coefficient, 

a = a + j (t^W -tW) log*. (2.25) 

Including both worldsheet instantons and the one-loop correction, the contact potential (equiva- 
lently the one-loop corrected four- dimensional dilaton) is now given by 



r = e * = ^ V - ^r^rl + e*« - (2.26) 

2 8(2tt) 3 2 192*' v ; 

2 

e*™ = -p— n K Re [ Li 3 (e 27TikaZa ) + 2nk a t a Li 2 ( e 2 ^ 2 ")] . (2.27) 

In addition to these perturbative corrections, one expects instanton corrections from Euclidean 
D(-l), Dl, D3, D5-branes wrapping supersymmetric 0,2,4,6-cycles on 2), as well as from Euclidean 
NS5-branes wrapping 2). Lacking a non-perturbative definition of string theory, an effective way 
of computing these effects is to assume that a certain discrete subgroup of the classical symmetry 
group SX(2,R) x N stays unbroken at the quantum level, and investigate how this symmetry 
constrains quantum contributions. 

For what concerns the nilpotent algebra of isometries N, it is natural to assume that quantum 
corrections will be invariant under the unipotent group iV(Z) consisting of integer monodromies 
around the large volume point, and large gauge transformations of the RR and NS axions. The 



2 



where 



Effective means k a G Z + for all a, not all fc 's vanishing simultaneously. 
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former consists of monodromies M e a in ( |2.9| ) with e a integer. The latter were analyzed in [Q, [JT 
using properties of the five-brane partition function. It was found that large gauge transformations 
act by integer shifts of the RR axions £ A , £a, but require an additional shift of the NS-axion, 



T, 



(r] A ,fjA,K) 



ct^o + 2k- 77 A (C A - 2d A ) + V A (( A - 



!>a) - V A Va 



(2.28) 



where 9 A , 0a are the characteristics entering the fivebrane partition function, and k is an arbitrary 
integer. Clearly, the same correction must also affect the transformation of the Darboux coordinate 
a in ( fl7D . 

Finally, since the SL(2, M.) action fl2.10D descends from the S-duality of ten-dimensional type IIB 
string theory (or, in the T-dual picture, from the diffeomorphism symmetry of M-theory compact- 
ified on 2) x T 2 ), it is natural to assume that the modular subgroup SL(2, Z) remains a symmetry 



at the quantum level. It was noted in ||31|| , and will be further confirmed in this paper, that the 
action of S-duality must however differ from the naive action ( |2. 10| ) by an additional shift in the 
RR axion dual to D3-branes 

c a ^c a - c 2 , a e(g) , (2.29) 
where e(g) G (— ~, ~] is defined in terms of the multiplier system of the Dedekind eta function i](t): 

2nie( a ) _ V (Sj) „=( ah \ 

' {CT + d) 1 ' 2 7]^)' y \cdj~ 

Of course, the same correction must affect the Darboux coordinate £' a , 



(2.30) 



C H- C + 



K a bcti b C ~ C 2 ,a e(flf) 



2(ce° + d) 



(2.31) 



These additional shifts are closely related to the quantization condition ( |2.7|) for the charge q' , and 
ensure the coincidence of the Heisenberg shift (° i— >■ (° + f]° with the SL(2, Z) transformation r i— > 
t + 77°, which was noticed at the classical level. Moreover, it was observed that the characteristics 
9 A must vanish, and that 4> a must transform under monodromies M e a as (fi a <p a + ^K aab e b . 

Based on these considerations, it is therefore natural to expect that the instanton corrected 
moduli space wil be invariant under the action of 5 SX(2,Z) k iV(Z), where the action of the 
unipotent group iV(Z) in the type IIB field basis is summarized in the following table: 





5b a 


5c a 




5 c 


5i/j 


T (l) 

(e a ,0) 









-e a c a - | n abc e a (b b + 2e b )c c 




T (l) 

(o,v a ) 





rf 


-\ K ab cV b b c + A abV b 


| K ahc r) a b b b c + A 0aV a 


V a (c a -(j) a ) + lA ab r} a ri b 
-lK abcV a b b (c c + 2r l c ) 


T (2) 
Va 








Va 








T (3) 
(vo,k) 











Vo 





In the remainder of this paper, we shall ignore the action of this discrete symmetry on the 
axions Co and ip, since they couple only to D5-brane and NS5-branes. In this sector, the generators 



'Further enhancements of S-duality to SX(3,Z) or SU(2, 1,Z) were proposed in 
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Tfa\ act trivially, and the remaining generators satisfy the Heisenberg-type relations (not to be 
confused with the Heisenberg algebra satisfied by the generators 



rpW rp(l) rp(l) rp(l) _ rp{2) qq\ 

1 (0,-v") ■ 1 (-e",0) ' 1 (0,^) ' 1 (e»,0) ~ 1 Kabc eV " K*' 6Z ) 

It is interesting to note that this Heisenberg-type algebra is isomorphic to the one which appears in 
the Coulomb branch of five-dimensional gauge theories compactified on a torus . We comment 
on the relation between the string theory and field theory set-ups in §|5|. 

2.5 Twistorial construction of D-instanton corrections 

Using the twistorial methods reviewed above and the symmetry restrictions discussed in the previ- 
ous subsection, in |27], ^] it was shown how to incorporate all D-brane multi-instanton corrections. 
For convenience we shall phrase this construction in the language of type IIA string theory compact- 
ified on the mirror threefold 2), where D-instantons correspond to Euclidean D2-branes wrapping 
special Lagrangian cycles with homology class 7 = (p A , q A ) G -£^3(2), Z) (or more generally, elements 
in the Fukaya category of 2J). 

To specify the twistor space corrected by these D2-brane contributions, one should replace the 
patch U , appearing in the above perturbative construction, by a (in general infinite) set of angular 
sectors which extend between two consecutive BPS rays 

£ 7 = {t g P 1 : Z^z^/t G HIT}, (2.33) 

where Zj(z a ) is the central charge 

Z 1 (z a ) = q A z A -p A F A (z a ), (2.34) 

evaluated at a fixed point z a in Msk- Across each BPS ray £ 7 , the Darboux coordinates are 
required to jump by the following complex contact transformation 

A£ A =^V log [1 - a^} 
Zm 

fi(V z a ) 

= V . ; QA log [1 - a^] (2.35) 

A « =^P- ( Li 2 K*7) - \ l0 g*7 [1 - ^7*7]) > 

where = E^o a |a — Qa^Ji ^(Ti z°) are the generalized Donaldson-Thomas (DT) invariants 6 , a 7 
is the quadratic refinement 

cr 7 = E(-~q A p A + q A 9 A - P a 0a) , (2.36) 



and we introduced the notation E(x) = e 2mx which will be extensively used below. In [28| it was 



shown how to integrate the contact transformation (|2.35| ) to a generating function = 



6 As further discussed in §3.1, the DT invariants are piecewise constant functions of the moduli z a . Thus, despite 
appearances ( [2.35 ) is still a holomorphic complex transformation. Moreover, as in the gauge theory setting pfj, the 
consistency of the construction for different values of the moduli z a requires that the jump of the DT invariants across 
lines of marginal stability (where two BPS rays collide) satisfies the Kontsevich-Soibelman wall-crossing formula Eq] . 
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While the result is somewhat cumbersome, we will only require its limit in the one-instanton 
approximation, 

H^,0 = ^^U 2 (a^). (2.37) 



Using the fact that the Darboux coordinates must approach the classical answer ( |2.15| ), ( |2.25| 



exponentially fast at t — 0, t = oo, the gluing conditions ( p.35| ) can be rewritten as integral equa- 
tions for the Fourier modes Af 7 , having the form of the Thermodynamical Bethe Ansatz equations 



MM 



Af 7 



< E E n (-r; *") (7, 7'} jT y ^ log [1 - ay Ay] j , (2.38) 



where 



A- = E(^fz 7 t-^)+A A -^ A ) (2.30) 



are the Fourier modes constructed from the perturbative Darboux coordinates ( 2.15|) . Eq. ( |2.38| ) 
uniquely determines the Darboux coordinates in each angular sector in terms of the DT invariants 
Q('j]z a ). In the weak coupling limit, it can be solved iteratively, leading to a (formal) multi- 
instanton series. The one-instanton approximation is obtained by substituting X 1 — > Ay on the 
r.h.s. and evaluating the integral over i! . The latter is dominated by a saddle point at 

* T = i r^T, S, = r 2 \Z,\ - i(/C A - g A C A ), (2.40) 

I 7 1 

leading to a correction of order 0(e~ 2nS ~>) consistent with the expected semi-classical action for 
single-centered D-instantons. Having solved these equations, the contact potential is then obtained 
from the Penrose-type integral 



e * = i e "'" > + iHf - els E a (r, z°) f j (r% - t~ Zl ) log [1 - . (2.41) 

7 "^v 



Since the HM moduli spaces of type IIA and type IIB string theories compactified on mirror 
Calabi-Yau threefolds must be identical, the above construction also defines in principle the D- 
instanton corrections in the type IIB formulation. In practice however, this is complicated by the 
fact that the relation between the coordinates z a ,C A ,C\,<7 appearing in the boundary conditions 
fl2.15|) and the physical type IIB fields 7 is a priori unknown beyond the classical mirror map (|2.6|) . 
One way to specify this relation is to require that the physical type IIB fields T,c a ,t a ,b a ,c a ,co,ilj 
should satisfy the same transformation rules (|2.10|) as at tree-level. But then, one faces the second 
shortcoming of the construction above, namely that of the lack of manifest S-duality invariance. 



This problem was solved in the Dl-D(-l) sector in [^5[ |36| where both, a manifestly S-duality 
invariant twistor construction and the quantum corrected mirror map were found. In this pa- 
per we address the more complicated case of D3-instantons and uncover their S-duality invariant 
description in the linearized approximation. 



7 In contrast, it is reasonable to assume that the coordinates z a , £ A , (a, o are identified with the physical type IIA 
fields, as they transform in the required way under symplectic transformations. 
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3. Modularity of DT invariants for D3-instantons 



In this section, we establish the modular properties satisfied by the generalized DT invariants 
which count D3-instantons. The arguments are very similar to those which allow to establish the 
modularity of the partition function of D4-D2-D0 black holes in type IIA/2J, or equivalently of 



M5-branes in M-theory on 2J |37|, [51], |52], ^3|. This is not a surprise, since T-duality relates the 
HM moduli space in type IIB on 2) x Si to the VM moduli space in type IIA on 2) x Si, and maps 
D3-brane instantons wrapping a divisor T> in type IIB to a D4-brane instanton wrapping T> x S 1 in 



type IIA. In §3.1| , we recall the geometric representation of D3-instantons as coherent sheaves on 
divisors, and recall the main properties of the generalized DT invariants which count them. This 
subsection may be skipped by the reader uninterested in mathematical aspects of D3-instantons. 



In § p.2| , we express the DT invariants in the large volume limit in terms of the MSW invariants, 
which have simple modularity properties, and define a novel multi-instanton expansion which is 
well suited for modularity. In § |3.3j , we give a heuristic description of D3-instanton corrections in 
twistor space, formally summing up the transition functions while ignoring the fact that each of 
them is associated to a different contour. A more rigorous treatment of modularity in twistor space 
is deferred till §[|. 

3.1 D3-instantons, coherent sheaves and DT invariants 

First, we recall the description of D3-brane instantons in type IIB string theory compactified on a 
Calabi-Yau threefold 2) in the large volume limit. Similarly to their D4-brane cousins [[54], |55|, [56], [)3|] , 
D3-instantons are described by pure, dimension-two, semi-stable sheaves £ supported on a divisor 
T> C 2), i.e. a complex hypersurface in 2) (more generally, objects in the derived category of 
coherent sheaves [^7], Q.). Let the 4-cycles 7 a (Poincare dual to 2-forms u a ) denote an integer 
basis of i?4(2),Z) = A, and the 2-cycles 7 a (Poincare dual to 4-forms u a ) denote an integer basis 
of F 2 (2), Z) = A*, such that 

u a A u b = K abc u c , co a Au b = 5 b a uj^ v / UJ b = uj a = 5 b 2 , (3.1) 

where is the volume form, normalized to = 1, and K abc = f^^a^b^c = (ja, 76, 7c)- 

The homology class of the divisor T> may be expanded on the basis of 4-cycles as T> = d a, j a , or 
equivalently its Poincare dual may be expanded on a basis of 2-forms, [V] = d a co a . We denote by d 
the vector (d 1 , . . . , d b2 ). We assume that T> is an ample divisor, i.e. that [D] belongs to the Kahler 
cone, 

d 3 = K abc d a d h d c > 0, r ■ d 2 = n abc d a d b r c > 0, q ■ d = q a d a > 0, (3.2) 

for all effective divisors r a ^ a G iJ4(2J,Z) and effective curves q a j a E i^2(2),Z). The intersection 
matrix of 2-cycles of an ample divisor D provides a natural quadratic form K abc d c on A, which has 
signature (1, b 2 — 1). This also provides the quadratic form K ab = (n abc d c )~ l of A*. In the following, 
we shall use this quadratic form to identify A as a sublattice of A*, k a k a = K abc p b k c . 

Given any coherent sheaf ^ on 2), the D-brane charges are given by the components of the 
generalized Mukai charge vector 7' on a basis of fp ven (2), Z), 



V = ch E sjTV® =p°+ p a oo a - q' a u a + </ Q u$, (3.3) 
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and satisfy the quantization conditions ( J2.7| ). D3-brane instantons correspond to dimension- two 
sheaves with p° = 0. Such a sheaf £ can be considered as the push forward £ = ^(J 7 ) of a coherent 
sheaf J 7 on T>, where i : T> — » 2) is the embedding of the surface X> in the Calabi-Yau. The Chern 
character ch £ is related to the Chern character of the sheaf J 7 on 2) by the Grothendieck-Riemann- 



Roch formula |58|, [59 



ch(£) = (ch(J-) Td(P)) / Td(2J) , (3.4) 

where is the push-forward in cohomology, obtained by first Poincare dualizing on £>, then pushing- 
forward in homology, then Poincare dualizing on 2) (thereby increasing the form degree by 2). 
The Chern character ch(J-") is given in terms of the rank r of J 7 and Chern classes c^J 7 ) by: 
ch(J 7 ) = r + ci(J 7 ) + ^ci(J 7 ) 2 — C2(J 7 ). We recall that the first Chern class of J 7 must satisfy 
^(J 7 ) G i/ 2 (2},Z) + ~Ci(A/"), where A/" is the normal bundle of T> in 2). The components of the 
charge vector 7' are then p| B(| 



P° = 0, p a = r ef , = - J t > a ) (d(^) + T - d(2?)) . 

fl& = jf^(ci(^ + ^)) a + ^X^)-rA(n 



(3.5) 



where t* denotes the pull-back map in cohomology, and ci(X>) = —c\(N) = —[D]. The Euler 
number of T> is 

X (V) = [Vf + c 2 (2J) ■ [V] = d 3 + c 2 • d, (3.6) 
while A (J 7 ) is the Bogomolov discriminant, invariant under tensoring by line bundles, 

1 f ( , ,^ J/ 



A(^) = - ^ ^ 2 (^) - cf(^) j . (3.7) 

Similarly to BPS D4-brane black holes, D3-instantons with the minimal number of fermionic 
zero modes (namely, four) correspond to semi-stable sheaves. Recall that a sheaf J 7 is semi-stable if 
all subsheaves J 7 ' C T satisfy ^(J 7 ') < ^(J 7 ). The (generalized) slope 8 for D3-instantons is defined 
by 

m= «±£l. (3.8) 
p ■ t 



In particular the Bogomolov discriminant ( p.7|) of a semi-stable sheaf is positive, which can be 
stated physically as a bound 



qo<^x(V) = ^(^ 2 p 3 + c 2 .p), (3.9) 



on the 'invariant D(-l) charge' 



^9i-^' 2 = ^(P)-rA(J), (3.10) 



8 The relevant stability condition in string theory is IT-stability, but this reduces in the large volume limit to 
(generalized) slope stability (see p3 for further details). 
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where q' 2 = K ab q' a q' b . A useful property of go is that it is, like the Bogomolov discriminant, invariant 
under tensoring of J 7 by a line bundle £ on V. If c\(C) = — e G A, the corresponding change of the 
Chern character ch(£) is the 'spectral flow' 

p^p, q'^q'-e, q' ^ % - e • q + - p ■ e 2 . (3.11) 

Note that the residue class of q' G A* modulo A is unaffected by this flow. This operation however 
does not leave invariant the stability condition for sheaves ^ on SI], unless one also changes the 
moduli z accordingly: z \-t z + e |38| . 

The moduli space of semi-stable dimension-two sheaves £ on 2) is a complex space of finite 
dimension. The generalized Donaldson-Thomas invariant Q(0,p,q',q' ; z) is an integer given, in- 
formally, by the Euler characteristics of this moduli space. It is also useful to consider the rational 



DT invariant [||, ||, |2j 



0(75 *) = X) ^ n frM *) • (3-12) 

which coincides with the integer / d; z) whenever 7 is a primitive vector, but in general is a 
rational number. Both Q and Cl are piecewise constant as a function of the Kahler moduli z, but 
are discontinuous across walls of marginal stability where the sheaf becomes unstable. This moduli 
dependence persists even in the large volume limit, with the exception of CY threefolds with 62 = 1, 
where the invariants fl(p, q' , q' ; z) are independent of z at large volume. 

3.2 MSW invariants and modularity 

The moduli dependence of the DT invariants makes it difficult to state their modular properties 
directly. Fortunately, it is possible to express the DT invariants in the large volume limit in terms 
of a different set of invariants Q p (q,q' ), which we call the MSW invariants. They are moduli 
independent and given by Fourier coefficients of a certain Jacobi form, namely the elliptic genus of 



the MSW superconformal field theory which describes D4-D2-D0 black holes |37|, tDJ I3> HI 
2 p (r, y) = Tr'(2J 3 ) 2 (-l) 2J3 e((l - |) r - (l - g) f + q' ■ y) 

= X (~ i r q dpW, Qo) e(-<?ot - \ q' 2 r - I q%f + q' ■ y) , (3.13) 

Here, the summation runs over 

q'GZ^ + ip, ?;6Z-ic r p, (3.14) 

where (p) a = n a bcP b p c , subject to the bound (|3.9|) . For a general vector k G A, the vectors 
k± G A ® R are projections of k onto the positive and negative definite subspaces of A <g) R defined 
by the magnetic charge vector p and the Kahler moduli 9 t, 

k ■ t 222 

kjf. = 7t k— k — kj^j k = k < H~ k (3.15) 

p ■ t 



9 This result is usually stated for t oc p, but the SCFT has a moduli space of marginal deformations which allows 
to go away from this point jO . 
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which satisfy 

k 2 + >0, k 2 _<0, Vfc^O. (3.16) 

We also use the notation k + to denote the modulus of the vector k + . It follows from general 
properties of the CFT that the elliptic genus ( |3.13| ) is a Jacobi form of weight (— |, ~) under 
SL(2, Z), with multiplier system Mz = E(e(g) c 2 • p), where e(g) is the multiplier system ( |2.30| ) of 
the Dedekind eta function. 

The MSW invariants Q p (q',q' ) coincide with the generalized DT invariants at the so-called 
'large volume attractor point' [Q 

ti p (q',q' )=ti(0,p,q',q' ;z oo ( 1 )), Zoo {j) = lim (b( 7 ) + iA t (7)) , (3.17) 

A— >+oo 

where 2(7) = 6(7) + it (7) is the standard attractor point. In particular, fi p (q',q' ) is invariant 
under tensoring by a line bundle ( |3.11| ), which corresponds to the spectral flow symmetry in the 
SCFT. Decomposing 

q' = /x + e + ±p, (3.18) 

where p, E A* /A is the residue class 10 of q' — |p modulo e, it follows that the MSW invariant 
&p(q' '? <?o) = &p,fj,{<lo) depends only on p, /x and %. As a result, the elliptic genus has a theta 
function decomposition: 

Z P (r,y,t)= ^2 h pA T )°pA T ^y^ t ) > ( 3 - 19 ) 

fj,eA*/A 

where 9 P)ll is the Siegel-Narain theta series 11 associated to the lattice A equipped with the quadratic 
form n a b of signature (1, 62 — 1), 

O w (T,y,t)= J2 (-l) kp v(^k 2 + T + ±k 2 _f + k- y y (3.20) 
keA+n+^p 

and the ^-independent coefficients are given by 

Vm(t)= E Ci p ^(q )E(-q r). (3.21) 

q Q <r X {T>)/2A 

Since the theta series ( p.20|) is a vector valued Jacobi form of modular weight (|, 3^) and multiplier 
system Mg under S 1,(2,7*), it follows that the function h Pll must transform as a vector- valued 
holomorphic 12 modular form of negative weight (— y — 1, 0) and multiplier system M(g) = Mz x 
Mg 1 under the full modular group SX(2,Z). This is equivalent to M(g) = Mz x since is 
unitary. The Fourier coefficients of h p M with g > are known as the polar degeneracies [53|, and 
determine the remaining MSW invariants by the usual Rademacher representation. 13 

10 The quotient A*/A is a finite group of order | detK a f,|, sometimes known as the discriminant group. 
11 Note that (3.20) coincides with the theta series ( |A.6| ) with y + = b + r — c + , y = fo_f — c , up to a phase 
V(\{b 2 + T + b 2 _?-b-c)). 

12 Non-compact directions in the target space of the CFT could potentially lead to mock modular forms and 
thus holomorphic anomalies pq |. For local Calabi-Yau manifolds, e.g. 0(—Kp-2) —> P 2 , it is known that the 
holomorphic generating series of DT-invariants for sheaves with rank > 1 requires a non-holomorphic addition in 
order to transform as a modular form |36]. We assume that this issue does not arise here. 

13 The most polar term q~ r x(' D )/ 24 j arises in the CFT from the fact that Lq > and that cl = r%(I>). Application 
of the Cardy formula shows that the asymptotic growth of states is dominated by r = 1 for given magnetic charge p, 
therefore the contribution of non-Abelian D3-instantons is exponentially suppressed compared to Abelian instantons. 
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Away from the large volume attractor point £00(7) (but still at large volume), the DT invariant 
0(0, p, q' , q' ; z) differs from the MSW invariant f2 PiM (go) by terms of higher order in the MSW 
invariants 



+ £ 5M^-W-w«tv»>) 

71+72 = (0,P,<j',?&) (3.22) 
Pi>0 V ' 

x ( 7l , 72 ) (_l)<™> Pl , Ml («Zo,i) ^W<7o, 2 ) 
+ . . . 

where 7« = (0, p i; q'j, g(, f ), and therefore (71, 72} = Pi ■ <? 2 — P2 ' 9i- Moreover one recognizes on the 
second line the difference of the slopes (|3.8|) of the constituents. The higher order terms can be 
thought of as describing bound states of the MSW constituents, which exist away from the large 
volume attractor point Zoo (7)- This decomposition is analogous to the decomposition of the index 
in terms of the multi-centered black hole bound states, although the elementary MSW constituents 
may themselves arise as bound states of more basic constituents such as D6-anti D6 bound states 
53j, which may decay across walls in the interior of the Kahler cone. The decomposition (|3.22| ) 



hinges on the fact that walls for constituents with p l = 0, p 2 7^ lie on the boundary or outside the 
Kahler cone, such that the only relevant walls are those between D3-instantons, i.e. constituents 
with charges p { > [p8 |. As shown in [38|, the 'two-centered' contribution in ( |3.22| ) (after smoothing 



out the sign function into a error function as in Eq. ( |A.8|) ) leads to a modular invariant partition 
function with the same modular properties as the elliptic genus, and it is expected that modularity 
persists to all orders. It is worthwhile to note that the expansion ( |3.22| ) has a finite number of 
terms in any chamber related by a finite number of wall-crossings from the large volume attractor 
chamber. 

Having expressed the generalized DT invariants in terms of the moduli-independent MSW 
invariants, which have simple modular properties, it is thus natural to reorganize the multi-instanton 
series expansion of the integral equations (|2.38|) by first expressing ft(j;z) in terms of Q Ptfl (q ) 



using ( 3.22| ), and then expanding in powers of the MSW invariants f2 PiM . In terms of the usual 



multi-instanton expansion of (|2.38|) , the n-th order in this reorganized multi-instanton expansion 
will contain the usual n-instanton contribution weighted by products of MSW degeneracies p/x , 
together with terms with n' <n instantons weighted by terms of order n— n! in the expansion ( |3.22|) . 
In the remainder of this work, we shall be concerned with only the first term in this expansion, 
which we dub the 'one-instanton approximation'. 

3.3 One-instanton approximation: heuristics 

In the one-instanton approximation introduced in the previous subsection, quantum corrections to 
the perturbative metric are effectively encoded by the set of holomorphic functions 14 

H^ , £, I) = J^y ^W<?0) CTyXy , (3.23) 



14 Indeed, the multi-covering contributions induced by the dilogarithm function in (2.37) are described just as well 



by replacing Li2(<r 7 A' 7 ) i— > a 1 X 1 and using the rational DT invariants ( 3.12 ) in place of the integer DT invariants. 
Due to (3.22), the rational DT invariants can then be replaced by O p , M . 
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which generate the infinitesimal complex contact transformation across the BPS ray £ 7 . The set of 
all (£y, H^) can be viewed as a Cech representative of H l (Z, 0(2)), which is known to classify the 
linear deformations of QK manifolds. In this subsection, we take a heuristic approach and consider 
the formal sum H = H^, ignoring the important fact that the functions H 1 are attached to 
different contours. This will allow us to highlight the main mechanism which ensures invariance 



under S-duality, deferring a rigorous treatment to § |4.1| onward. 

Ignoring the contributions of D5 and NS5-instantons, which are exponentially suppressed at 
large volume with respect to D3-Dl-D(-l)-instantons, and restricting to a fixed D3-brane charge 
p, the D3-D1-D(-1) one-instanton corrections are thus described by the formal sum 15 

i 

H V = T^A2 E *W?o) E(p l'-q' t- q' o e) , (3-24) 



where £\> iM (go) are the MSW invariants introduced in the previous subsection. In writing (|3.24j ), we 
identified the quadratic refinement cr 7 with the phase (— l) p q ' by setting the characteristics 9 A , a 
to zero. We also omit the phase E(|Aae£> A £> e ) which is unimportant for fixed magnetic charges 
p A . The assumption 0a = is for simplicity, since 0a couples only to magnetic charges which are 
kept constant in our set-up. In contrast, the condition 9 A = is necessary for modular invariance, 
which requires the phase factor (— l) p q ' in the sum ( p. 24 ). The vanishing of 6 A was also seen to be 



necessary in the linear analysis of NS5-instanton corrections in |31 

Using the same steps as before, we can therefore trade the sum over q' , q' in ( |3.24j ) for a sum 
over /x, % and fc = e + /x + |p, and represent ( |3.24| ) as a sum of theta series analogous to Q3.19| ), 

^> = 7^ E M') E ^(£°)H PiM (£ ,£), (3-25) 



(2tt) 
where 



Sp !M (£ ,£)= E ^{A 1 ^-V 2 -k-i-\ee\ j (3.26) 

and /i P)At (£°) is the same modular function which was defined in ( |3.21| ), evaluated at r = £°. It 
therefore transforms as a vector- valued holomorphic modular form of modular weight — y — 1 
and multiplier system M(g) = Mz x Mg under the S-duality action (2.2U|). On the other hand, 



Sp. /X (£ ,^) is formally a holomorphic theta series for the lattice A of rank 62 with quadratic form 
—K a b, with modular parameter r = £° and elliptic parameters £ a . Thus, under the action (|2.20|) , 



it is expected to transform as a holomorphic Jacobi form of weight y, multiplier system Mg 1 
and index m a b = —^k^. Finally, thanks to the term proportional to C2 ja in the transformation 
of the exponential prefactor in ( |3.25| ) transforms as the automorphy factor of a multi- variable 
holomorphic Jacobi theta series with the index m a fe = \n ab and multiplier system M% . We thus 
conclude that, under the action of S-duality in twistor space, the formal sum ( p.24j) transforms as 
a holomorphic Jacobi form of weight —1 and trivial multiplier system. In [[JBJ it was shown that 
this is precisely the condition ensuring the presence of the SL(2, Z) isometry group. Thus, H m is 
expected to generate an S-duality invariant deformation of the perturbative HM metric. 



15 



The symbol ^2 indicates a divergent sum. 
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However, this analysis overlooks the important fact that the quadratic form n ab has indefinite 
signature (1,62 — 1), and therefore the theta series ( |3.26| ) is formally divergent. Fortunately, the 
theta series ( |3.26| ) never arises as such in the computation of the metric, rather each of the terms in 
( |3.26| ) must be integrated along a different contour, which renders the resulting series convergent. 
Another fortunate circumstance is that there does exist a natural holomorphic 'mock' theta series 
of signature (1,&2 — 1), the convergence of which is ensured by restricting the sum to lattice 
vectors lying in a certain cone where the quadratic form is definite positive [KJ. While each 
of these two series transform with modular anomalies, their sum is in fact modular invariant, 
ensuring the modular invariance of the D3-instanton corrected metric, as the next section will aim 
to demonstrate. 



4. D3-instantons, period integrals and mock theta series 

In this section, we discuss how S-duality is realized in the presence of D3-D1-D(-1) instantons, but 
in the absence of D5 and NS5-instantons. This situation arises in the large volume limit, where 
D5 and NS5-brane instantons are suppressed compared to D3-D1-D(-1) and worldsheet instantons. 



In § |4.1| we introduce a large volume limit which retains modular invariance but reduces the D3- 
instanton sums to (generalized) Gaussian theta series, and obtain the instanton corrections to the 
Darboux coordinates and contact potential in this limit. We establish their modularity (or lack 



thereof) in § Q , §4.3| . Finally, in § |4.4| we show how modular anomalies can be canceled by a contact 



transformation generated by a mock theta series. 

4.1 Instanton corrections in the large volume limit 

In order to analyze the action of S-duality on the D-instanton corrected twistor space, it is useful to 
focus on the Darboux coordinates around the point t = — i, or equivalently z — 0, which is invariant 
under the S-duality action ( |2.21| ) on the P 1 fiber. In the large volume limit t a — > 00, it may be 



checked that the saddle points ( |2.40| ) dominating the integral equations ( |2.38| ), for D3-D1-D(-1) 
charges 7 = (0,p, q, go) such that p ■ t 2 > 0, accumulate near the same point 

«>-.-2<*+% + ..., <*-i^±% + .... (4,, 
vP ■ t VP- * 

Thus, it is natural to expand the Darboux coordinates around the point z = 0, t a — > 00, keeping 



the product zt a fixed. Retaining all the terms from the perturbative result ( |2.15| ) which are finite 
or divergent in this limit, we find 

£ a = rb a - c a + 2r 2 t a z + 5£ a , 

i' a = l -r 2 K abc t b t c + ~c a + -n abc b\c c - rb c ) - 2r 2 n abc t h b c z + ir 2 n abc t b t c z 2 + 5% a , 

Co = \^abcT2t a t\2t C Z + 3b c ) +5 - \n abc b a b\ C C - Tb C ) + T 2 K^b^b* Z ~ itV) + 

o 

ol = -l^ abc r 2 t a t\2t c zf + 3c c ) + ^ + l^ ab cb a (b b T - c b )(b c r - 2c c ) (4.2) 
o o 



+T 2 zn abc t a (b b {b c r - 2c c ) + t b {2zb c r 2 + \z 2 c c )) + 5a', 



-20- 



where 5£ a , etc, denote the instanton corrections which include two types of contributions: from the 
integral terms in ( |2.38| ) and from instanton corrections to the classical mirror map ( |2.6| ). Note that 



the former contributions can be split into two parts. The first corresponds to D-branes with charge 
7 with p ■ t 2 > 0, while the second comes from anti D-branes with opposite charge —7. In the 
former case, the integrals are dominated by a saddle point at t' , which approaches if = — i in the 
large volume limit, while in the latter case, they are dominated by a saddle point at the antipodal 
point if = —1/t' which approaches t' = i or z' = 00 in the large volume limit. This suggests to 
keep fixed inside these integrals a different combination, t a /z'. 

Besides the Darboux coordinates, it is important to consider also the contact potential e*. In 
our approximation, it is dominated by the classical term in (|2.26|), so we can write it as 



= \r 2 2 V + 5e\ (4.3) 

Our task will be to compute the corrections <5£ a , etc. in the one-instanton and large volume 
approximation, and to show that the resulting instanton corrected Darboux coordinates (|4.2j ) and 
contact potential ( |4.3|) continue to transform in the same way ( |2.20|) , ( [2.23|) as their perturbative 



counterpart, possibly up to a contact transformation. For what concerns the coordinates £ a and 
the contact potential e*, this requires showing 

XCa $ 

6? H- 5e* h. j—- (4.4) 

ct + a \ct + a\ 

For what concerns the other coordinates £' a , £' , a', the corresponding constraints are best expressed 
by forming suitable linear combinations 16 

S£'a=S£' a + Kabc {b b -izt b ) 5?, 

5+a' = 5a' + r 8& + n abc (b a - izt a )(c b - rb b - 2T 2 zt a )5£ c , (4.5) 
5_a' = 5a' + f5£ + n abc (b a - izt a )(c b - nb b - r 2 zt a )5C- 

These linear combinations are designed to transform canonically under S-duality as 

Xt' v St' X ' v X ' v ^- a ' (A R\ 

ct + a ct + a 



if and only if the corrected Darboux coordinates ( |4.2|) transform according to ([2.20 ). 



As was mentioned above, the instanton corrections acquire two types of contributions. The 
first is computed by replacing X 1 by on the r.h.s. of the integral equations ( |2.38| ) and ( |2.41| ), 



taking the large volume limit, and evaluating the integrals in terms of the type IIA coordinates 
r 2 , z a , ( A , (a, <J- The latter can then be replaced by the type IIB coordinates T,t a ,b a ,c a ,c a ,co,ip, 
adapted for S-duality, using the classical mirror map ( |2.6|) . However, there are further contributions 
stemming from instanton corrections to this mirror map and coming from the tree level terms on 
the type IIA side. In Appendix [B] we extend the procedure of [36] and compute these corrections 



to the mirror map in our approximation. The result is presented explicitly in ( p,15|) . In the same 



16 



These linear combinations are equal to the large volume limit of the linear combinations in Eq. (B.6) of pq] , up 



to a divergent term which also transforms as in (4.6) and can therefore be dropped 
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appendix it is shown that, adding up the two types of contributions, the instanton corrections to 
the Darboux coordinates and the contact potential can be written as 

5£ = 2-iripJp, 
-D ■ J p , 



5 + a' 
8-a! 

6e* 



To 

—D D — 2 

L 7T 



(4.7) 



r a J 

[47rr 2 D_i + (c - fb) -D]-J p + ir 2 z(p ■ t 2 ) ^ / — - H 



-7' 



9A 



g + - (k + b - iz't) ■ (k + b - 3iz't) 



H~ + c.c. 



where 



2vri(^ - z) 7 ' 



(4- 



and P-i, -D, -D are the covariant derivative operators 



17 



1 



27Ti 



7 6»r 



(1 



2ir 2 



D = -—[d h + fd c -m{c-fb)} 
2t 2 



D = (D) 



(4.9) 



Moreover, H± 1 is the large volume limit of ( |3.23|) where one keeps fixed z ±l t a , respectively. Ex- 
plicitly, these two functions are given by 



H. 



H 



(_l) P k 



n p ^(q ) E iS cl - -(k + b) 2 + f - (q + -(k + b) 2 _ r + c ■ (fc + -6) + iQ 7 (*' 



7 



(27T) 



fi P)M (g ) E i£ci + ~(fc + b)%r + Uo + + b)t f - c • (k + -b) + \Q 



(4.10) 



where S c \ is the leading part of the Euclidean D3-brane action in the large volume limit, and Q±-y(z) 
is the only part which depends on z, 



r-2 



Sci = — P ■ t - ic ■ p, Q ±7 (z') =r 2 p-t \(z) ±i 



„ ±1 , ; (k + b)+ 



(4.11) 



The notation j £ signifies that each term must be integrated on a contour which extends from 



z' = — oo to z' = +oo and goes through the saddle point z' = —i(k + b) + / y p -t . This corresponds 
to the BPS ray £ 7 in the large volume limit. 

As we shall see momentarily, the corrections fl4.7|) do not quite transform as required in ( |4.6| ), 
but the modular anomaly can be absorbed by a contact transformation. This issue however does 
not arise in the case of the contact potential, which we consider first. 



17 Acting on theta functions of weight ((), fj), the operators 2?f,, D,D raise the modular weight by (2,0), (1,0), (0,1), 
respectively. To avoid cluttering, we abuse notation and declare that T>^ annihilates the classical, modular invariant 
contribution S c \ in i? 7 . 
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4.2 Modular invariance of the contact potential 



The one-instanton correction to the contact potential in ( f4.7| ) involves a Gaussian integral over z' 
(renamed as z) which is easily calculated. Rewriting the integrand 



5e 



dz 



q + -{k + b): 



1 



and dropping the total derivative, we arrive at 
$ _ t 2 



71 T 2 47TT2 



d z z + 



1 



lllT 2 P ■ t 



2 "z 



5e 



16n^2r 2 p-t 2 Qa 



g + ~(k + b) 2 _--^- 

2 8/TT 2 



x E(iS cl -~(k + b) 2 + f- (q + l(k + b)A r + c-(k + h) 



H^ + c.c. (4.12) 



c.c. 



(4.13) 



The sum over q\ is recognized as the partition function ( |3.19| ), except for the insertion in the square 
bracket. The latter in turn arises by acting with the covariant derivative T)_a, obtaining 



6e* 



t~2 e 



-2tt5 c1 



_V_3 h p ^(T)6 p ^(T,t,b,c) + c.c. 



(4.14) 



3 1' 
'2' 2> 



to 



I I) 

2' 2>i 



while 



The action of T>_3 on /i P)At ('r)6 l p iAt (r, i, b, c) raises the modular weight from 
the overall factor of r 2 reduces this to (— |, — |) (recall that the combination r 2 p ■ t 2 is modular 
invariant). Therefore, the correction transforms as required in ( |4.4|) . 



4.3 Instanton corrections to Darboux coordinates and period integrals 

We now turn to the instanton corrections to the Darboux coordinates. Unlike the contact potential, 
which is modular covariant, it will turn out that the Darboux coordinates have anomalous modular 
transformations, due to the fact that the contour £ 7 is not invariant. However, we shall be able 
to give a precise characterization of this modular anomaly, by rewriting the Penrose-type integral 
over z' as an Eichler integral. 

For simplicity we start from the last term in 5-a' ( |4.7| ). It is given by a Gaussian integral in 
1/z' which can be easily evaluated. The result is 



/ 2\ f dz' TT iz 

ir 2 z(p • O 2^ / j^ H ^ = 



71^ 



^2r 2 p-t 2 e- 2 ^ 



1 E V/*' 

fj,eA*/A 



T)6 Ptlt (T,t,b,c). (4.15) 



It transforms as a modular function of weight (0,-1) in full agreement with the required transfor- 
mation properties (|4.6|) . 

All the other instanton corrections to Darboux coordinates in the large volume limit are de- 
termined by the function J p (z) ( f4.8| ). Therefore, our prime interest is to understand its modular 
properties. To this aim, let us consider the integral appearing in this function, namely 



X{z) 



dz' 
z' — z 



E(iQ 7 (;0). 



(4.16) 
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By shifting the contour, this may be rewritten as 

dz' 



I(z) 



z — z — ict 



a 



(k + b). 



, P = ^27ir 2 p-t 2 . 



(4.17) 



Since z can be absorbed by a shift of a, we shall first consider the value of X at z = 0. Using the 
identity 



dz' 



z — \a 



-/3 2 z 2 



ivr sgn( Re (a)) e a2 ^ Erfc ( sgn( Re (a/3)) a/3 



valid for a, (3 G C with Re ((3 2 ) > and Re (a) ^ 0, the integral representation 

POO 

Erfc( v / ^a;) = / du u~ 1/2 e~™, Re (x) > 0, 



(4.18) 



(4.19) 



and changing variable from u to w = t — iu / (k + b) 2 +) valued in the lower half plane, one may cast 
X(0) in the form of an Eichler integral (|A.14 ) 18 



1(0) = -n 



dw 



W — T) 



(k + b) + E -(k + b) 2 + (f - w] 



(4.21) 



Inserting the representation ( 4.21|) into (4.8), restricting to z = for simplicity, and carrying 
out the sum over q\ in the same fashion as in §|3.2|, we arrive at 



l e 



-2nS c i 



/xeA*/A 



dw 



'i{w — T) 

x (-l)^(fc + fe) + E(-|(fc + 6)^-i(fc + 6) 2 r + c-(fc + i6)). 

fceA+/x+ip 



(4.22) 



The second line in this expression is recognized as the (complex conjugate of the) Siegel-Narain 
theta series ( |A.16|) of weight (|, \ (b 2 — 1)), analytically continued away form the slice w — f. Eq. 
( |4.22| ) is thus an Eichler integral of weight (y,0), multiplied by the vector- valued modular form 
h p ^ of weight (—1 — y, 0). Thus, J v (z = 0) transforms with modular weight (—1, 0), except for a 
modular anomaly of the form ( A.15|) given by the period integral of the theta series ( |A.16| ). 

More generally, using this type of manipulations, one may rewrite the full function J v (z) away 
from z = as 

j g— 27rS c i _^ 



peA*/A 



T M (w,r;2:) dw 



(4.23) 



W — T) 



In fact, the integral over z' (4.18) can be rewritten as an integral over w without ever evaluating it in terms of 
the error function, upon representing both factors in the integrand in terms of their Fourier transforms 

1 



i sgn(a) 



z' — \a 



7?(sgn(a) u) 



(4.20) 



" 2|/3|V5F. 

where r\ is the Heaviside step function, and carrying out the integral over z 1 ,uj,uj' . Thus, the integration variable w 
in (4.21) is Fourier dual to the twistor coordinate z' in (4.16). 
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where 



T tl (w,f;z)=E^r 2 p-t 2 ^—^z 2 ^ £ (-I)** Uk + b) + + 



fe6 A" + i P V W-T J (424) 

x E(-|(fe + bf + w - \{k + 6) 2 r + c-(k + ±6)) . 



It is easy to see that the integral is identical to the complex conjugate of the generating function 



&p lt j,(y) ( IA..27D evaluated at y = 2r 2 ^/p • t 2 z, i.e. 



J M = -MT E KA^QtA^Vp^z). (4.25) 

/x6A*/A 



In Appendix ^ it is shown that provided y transforms with the modular weight (0, —1), which is 
indeed the case for our identification between y and z, the generating function QpJ^y) transforms 
as an Eichler integral of weight (0, 6 2 /2). As a result, the full function Q4.25|) possesses the following 
transformation property 



JM - (- + d)- [ JM,) - — £ ^(r) /_ ^ j , (4.26) 



i.e. transforms with modular weight (—1,0) and an anomaly given by a period integral. Recalling 
the properties of the covariant derivative operators Pf,, D, D mentioned in footnote this allows 
to conclude that the corrections 5£, 5£, 8±a' transform as required in ( |4.6|) , except for the modular 
anomaly ([4.26 ), acted upon by the corresponding covariant derivatives. 



4.4 Anomaly cancelation from mock theta series 

Having expressed the corrections to the Darboux coordinates as Eichler integrals, or modular 
derivative thereof, it remains to show how the modular anomalies can be canceled by a suitable 
contact transformation. 

As explained in Appendix the modular anomaly of an Eichler integral can be canceled by 
adding to it a certain mock theta series. Typically it appears as a theta series with insertion of a 
difference of two sign functions defined by two Kahler moduli t and t', with one of them chosen to 
lie on the boundary of the Kahler cone, i.e. pi! 2 = 0. The simplest example is provided by ( |A.ll ) 
and ( |A.10| ), respectively, which describe the instanton corrections to J7p(0). 

Similarly, the modular anomaly in (|4.26|) can be canceled by replacing the generating function 
Qp in ( f4.25|) by its modular completion — £?® )M , where the second term is the generating 



function defined in (|A.24j ), evaluated at y — 2r 2 yp- t 2 z. In this way, we find that the modular 
completion of J p {z) is given by 

Cfp iJp ^anom > (4.27) 



-2ttS c1 



where 

^o m = ^ r h p ^r)g^(2r 2 ^/p~¥z) . (4.28) 

neA*/A 
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Using ( A.26| ), this can be rewritten as 



tieA*/A 

("l) fc - p [sgn ((fc + 6) • t) - sgn ((fc + 6) • f )] E (p.£-fc.£-^° (*0 2 ) , 



x 



which shows that miraculously, is a holomorphic function of the Darboux coordinates. In 

fact, H anom is proportional to the holomorphic mock theta series QA.10 ) evaluated at r = £°, y = £ 



(upon taking into account the identification between y and 6, c mentioned in footnote |TT|). 

Moreover, one can show that the corresponding modular completions of the Darboux coordi- 
nates 5£,5£' ,5a', following from (|4.5|) , (|4.7|) and obtained by substitution of J v by J7" p , have the 



following form 

<^-^#ano m , C + fe' + (1 - £ A %) # anom (4.30) 

and transform according to ( |2.20|) , without any further anomaly. Moreover, the modular completion 
( 4.30Q is recognized as a holomorphic contact transformation ( |2.13|) . It is interesting to note that 



the theta series Q4.29|) can be written as 

#anom = ~ J] [sgn {{k + b) ■ t) - Sgn {{k + b) • *')] #7 ( 4 '31) 



2 



and differs from the formal sum ( |3.24|) only by an insertion of the two signs which makes the sum over 



charges convergent. As a result, the compensating contact transformation is generated by a proper 
subset of the original transition functions, and the modular covariant Darboux coordinates differ 
from the original type IIA Darboux coordinates by a sequence of symplectomorphisms associated 
to the BPS states for which the two signs are different. This mimicks the situation in the Dl-D(-l) 



sector 36 



Thus, we have demonstrated that the modular anomaly is canceled by performing the trans- 
formation ( p,13|) generated by the mock theta series (|4.29| ). In fact, we have a continuous family 
of modular completions labeled by the parameter t' on the boundary of the Kahler cone. It would 
be interesting if this ambiguity was fixed by some physical principle. This concludes the proof of 
the modular invariance of the HM moduli space corrected by D3-instantons in the large volume, 
one-instanton approximation. 



5. Discussion 



In this work, we investigated the modular invariance of D3-Dl-D(-l)-instanton corrections to the 
hypermultiplet moduli space of type IIB string theory compactified on a CY threefold 2). We 
focused on the large volume limit, where D5-NS5-instantons can be consistently ignored. Using 
similar arguments as the ones which enter in the proof of modular invariance of the partition 
function of D4-D2-D0 black holes in type IIA string theory compactified on the same threefold 2), 
we showed that the DT invariants which govern D3-instanton corrections are expressible in terms 
of the Fourier coefficients of the elliptic genus of the MSW superconformal field theory, which we 
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refer to as MSW invariants. Unlike the DT invariants, the MSW invariants are independent of 
the moduli and have simple modular properties. In the one-instanton approximation, we found 
that the D3-instanton corrections to the standard Darboux coordinates are expressible as Eichler 
integrals of the MSW elliptic genus, and modular derivatives thereof. Thus, they suffer from 
modular anomalies, which can however be absorbed by a local complex contact transformation. 
The generating function for this contact transformation was recognized as the holomorphic mock 



theta series of signature (1,62 — 1) introduced in [40]. In this physical set-up, the Eichler integral 



which provides the modular completion of the mock theta series arises as a Penrose-type integral 
along the twistor fiber, or rather its Fourier transform. 19 Unlike Darboux coordinates, the contact 
potential (also known as the four- dimensional dilaton) transforms covariantly, and is proportional 
to the (modular derivative of the) MSW elliptic genus. 

While it is very satisfying to see modularity emerge in the large volume, one-instanton ap- 
proximation, an outstanding problem is to extend the results of this paper beyond this regime. 
This can be in principle addressed by expanding the integral equations (|2.38| ) to higher order in 



the MSW invariants, as indicated in §3.2j . There are two contributions at next-to-leading order, 
namely (i) the two-instanton correction in the usual iterative expansion of ( |2.38|) , with the DT 
invariant ^(7; z) replaced by the MSW invariant f2 PiM , and (ii) the usual one-instanton correc- 
tion, with 0(7; z) replaced by the quadratic term in ( |3.22| ). The first contribution is weighted by 
e -2wr 2 (\z( 7l ,z)\+\z( 72 ,z)\) ^ w hereas the second contribution is weighted by e - 27rT2 l z ^ 1+ ^' z )l. By con- 
struction, the sum of these contributions is continuous across walls of marginal stability, although 
they are separately discontinuous. We expect that the D3-instanton correction at two-instanton 
order will be modular invariant, and will be related to the two-center D4-brane partition function in 
the same fashion as at one-instanton order, with the two-center Siegel-Narain theta series replaced 
by an indefinite theta series. It is encouraging that the partition function for two-center D4-black 
holes, corresponding to the contribution of (i), indeed becomes modular invariant and continuous 
after adding a suitable non-holomorphic completion ||38|| , which should arise from contributions of 
type (ii). 

More ambitiously, it would be very interesting to show that S-duality invariance holds at the 
non-linear level. One way to achieve this would be to recast the standard type IIA twistorial 



construction into the manifestly invariant framework developed in [[36]. A naive attempt using 
the standard representation of the MSW elliptic genus as a Poincare series [42| fails, due to the 
necessity of regulating the Poincare series (see Appendix y). It is conceivable that this obstruction 
may be avoided by taking into account the constraints on the polar degeneracies from modular 
invariance. Eventually, we hope that it will be possible to formulate all D5-D3-D1-D(-1) and NS5- 
brane instanton corrections in a manifestly invariant S-duality fashion, going beyond the linear 



order analysis of |3T[] and providing a tight set of constraints on the generalized Donaldson-Thomas 
invariants on Calabi-Yau threefolds. 

Finally, we expect that the structure uncovered in this paper in the context of the HM moduli 
space in type II string theories on a compact CY threefold will continue to exist in the rigid limit 
near singularities of 2). Specifically, the Coulomb branch of M = 2 five- dimensional SU(2) gauge 



theories compactified to a two-torus p5 |, obtained by considering type IIA string theory on a non- 



19 See [p5[ for another example where the modular completion of a mock modular form arises naturally from 
physics. 
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compact CY threefold given by the anti-canonical bundle O(-Kg) — > S over a rational surface S 
[ 57] , gives a family of hyperkahler manifolds parametrized by the modular parameter of the torus, 
which must be invariant under modular transformations. We expect this family of HK metrics 
to be described by a similar twistorial construction as in this paper, using the generalization of 
the QK/HK correspondence put forward in ||36|| . For such line bundles over rational surfaces, 
the DT invariants 0(7; z a ) with p° = can be computed explicitly f68 |, and the analogue of the 
MSW CFT is a (0,4)-supersymmetric sigma model with target space given by the Atiyah-Hitchin 
moduli space of monopoles in three dimensions p5| j9[| . An important difference between the rigid 



case and the assumptions in this paper is that the generating function hp^ij) of DT invariants 
for rational surfaces is mock modular [^6], 70|. It is natural to expect that its non-holomorphic 
modular completion will arise from multi-instanton corrections. Finally, by the general logic of 
the QK/HK correspondence, this family of HK metrics will be naturally endowed with a rank two 
modular invariant hyperholomorphic line bundle analogous to the rank one hyperholomorphic line 



bundle which arises in 4D gauge theories fl30| , 71 1, whose physical role remains to be understood. 



Acknowledgments 



We would like to thank Daniel Persson for valuable comments on the draft. 



A. Indefinite theta series and period integrals 

In this Appendix, we review some aspects of indefinite theta series with Lorentzian signature 
— 1), which are essential for describing D3-instanton corrections. 

A.l Vigneras' theorem 

Let Q{x) = x 2 be a quadratic form on M. n with signature (1,ti), and A an n-dimensional lattice 
such that Q(k) G Z for k G A. Let V[x) be a function on M. n such that V(x) e nx2 is integrable. Let 
p G A be a characteristic vector (such that k 2 + k ■ p G 2Z, VkGA),^GA*/Aa glue vector, and 
A an arbitrary integer. Following [72] , with some changes of notations, we construct the family of 
theta series 

^ (V, A; r; b, c) = r" A/2 £ {-l) h *> V {^/W 2 {k + b)) E(-c ■ (k + ±6)) , (A.l) 



where q = e 2nlT and b, c 6 A £g> R. Irrespective of the choice of V and A, ([A.l ) satisfies the elliptic 
transformation properties 

P ,„ (V, A; r; b + k, c) =(-l) k ^ E(|c ■ k) ^ (V, A; r; b, c) , 
$ P ,„ (V, A; r;b,c + k) =(-l) k '*> E(-|& • k) ^ (V, A; r; 6, c) , 

for any k G A. 
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Now, let us assume that V(x) e nx2 decays sufficiently fast 20 at infinity, and that V(x) is anni- 
hilated by the differential operator 

[A + 2n(x-d x - \)]V(x) = 0, (A.3) 



where A = cr is the Laplacian operator on W 1 . Then, using the same methods as in f72|, one may 



show that the family of theta series ( |A.1| ), with running over all cosets in A* /A, transforms as a 
vector-valued Jacobi form of weight (0, A + n/2) under the full modular group. This means that it 
satisfies the modular properties 

^(P,A;-l/r;c,-6) = ^==(if) A+ fEf-lp 2 ') £ E(-/i • u) tf PjM (V, A; r; b, c) , 

J A* /A V 4 / ,JT7 /A 



i/SA*/A 



^ p ,m (P, A; r + 1; 6, c + 6)) = E(|( M + \pf) ^ (V, A; r; b, c) , (A.4) 

on top of the elliptic properties ( |A.2|) . In the next subsections, we present the relevant instances 
of this general construction relevant for the present work. 

A. 2 Siegel-Narain Theta series 

The simplest example arises for V{x) = e~^ x + 1 , A = -1, where x + is the projection of x along 
a fixed time-like vector t in W 1 (t 2 > 0): 

xf = ^±, x® = ^t = x-x® (A.5) 
yt t t • t 

so that (x^) 2 = (x^) 2 > 0, (x^) 2 < 0. In case no confusion arises, we write x+} = x + . The 
corresponding theta series ( |A.1| ) is then, up to a power of T2 and lattice shift, the standard Siegel- 
Narain theta series, 



9 Pjfl (r, t, b, c) =t~ 1/2 ^ ( e"™+, I : r: b. c 



(-l) k p q2 {k+b) + q-2 {k+b) - E(-c ■ (k + §6)) . 



(A.6) 



Thus, 6 P41 transforms as a vector- valued Jacobi form of weight (|, 
A.3 Zwegers' mock theta series 

The other prime example is Zwegers' indefinite theta series of weight (0, ~) |4(| (based on earlier 
works 0,0), 

p , M (r, t, t', b, c) = ^ [Erf (V2^(k + 6)?) - Erf (V2^(k + 

fceA+ M +|p (A.7) 

xg-^ 2 E(-c.(fc + |)), 



20 More precisely, f(x) = V(x) e* x should be such that, for any D(x) be any differential operator of order < 2 and 
R(x) any polynomial of degree < 2, f(x), D(x)f(x) and R(x)f(x) should be both integrable and square integrable 
on R n . 
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where t and t' are two time-like vectors and Erf (y/nx) = 2 f*e nt2 dt is the standard error function. 
This series again belongs to the class flOD with V{x) = Erf( v ^ra;J ) ) - Erf( v ^ra;J' ) ). Each of the 
two terms in V(x) satisfies the differential equation (|A.3|) for A = 0, and the sum of the two terms 
satisfies the decay conditions, although this is not true for each term separately. Using 

Erf(v^Fx) = sgn(x) (l - Erfc(v^F|x|)) , (A.8) 

one may decompose (|A.7|) into a sum of three terms 

6 Pi/i (r, t, t', b, c) = -$ PlM (r, t, b, c) + 9 p , m (t, t, t', b, c) + $ p>M (r, t', b, c), (A.9) 

where the middle term is holomorphic in r, 



E 



vfep 



sgn((A; + b) { ?) - sgn((A; + b)^) q-^ k+b ^ E(-c ■ (fc + §&)) , 



»(*0 



while the first and last term are given by 



(A.10) 



^ (-l) fc - p sgn((A;+6)J ) ) Erfc (v 7 ^!^ + g-^ (fc+b)2 E(-c • (k + §&)) 



fceA+p+gP 



(A.ll) 

It is useful to note that both ( A. 10 ) and ( A.ll ) are of the form ( A.l ) with a function V(x) which is 
square integrable 21 and obeys the differential equation (|A.3|) everywhere, except on the locus where 
x+ vanishes. As a result, the three terms in ( |A.9|) are not separately modular invariant, although 
their sum is. 

The modular anomaly of ( |A.11|) may be exposed by using the identity 

(k + b) ( * ] dw 

yj-l{w ~ f) 

(A.12) 



sgn((A; + b)f) e 2 ^ k+h) + Erfc (V2^\(k + 6)J } |) = -i jf°° E^(fe + b)\{w - r)J 



to rewrite ( |A.11| ) clS db (generalized) Eichler integral, 

•V,(r,t,M = -ljT E _ (-■)% iWft P W1 + 



(A.13) 

where qv, = E(u>). Recall that if F(r, f ) is an analytic modular form of weight (f), f)), and if F(w, f) 
is its analytic continuation away from the slice w — r, the Eichler integral 



F(w, r) dw 



(A. 14) 



transforms with modular weight (0, f) + 2 — f)), up to modular anomaly given by a period integral, 

3 F(w,f)dw 



$(7r) = (cr + d) 



f)+2-f) 



$(r) 



-d/c 



[-i(w - r)] 2 -l 



(A.15) 



21 For the holomorphic theta series ( A.lCf ), this follows from the fact that the difference of sign functions vanishes 
on all time-like vectors. 
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The function F(r,f) is known as the shadow of the Eichler integral ( A.14Q , and can be extracted 
from $(t) by acting with the operator — (2t2) 2_ ^<9 t . 

Returning to ( |A.13| ), $ PjM (t, t, b, c) is then identified as an Eichler integral of weight (0,1), 
with shadow proportional to the theta series of weight (|, ^^r), 



r. 



2) e" 



-, -2;r; b, c 



^ (A; + 6) + g2( fc+b )+ g"2( fe + 6 )- E(-C ■ (Jfe + §&)) . 



(A.16) 



fceA+^+^p 



Let us now consider what happens with ( |A.9| ) when one of the time-like vectors, say ap- 
proaches the boundary of the Kahler cone, i.e. becomes light-like t' 2 = 0. Then (k + by* diverges 
and, since Erfc(|x|) is exponentially suppressed as x — > oo, the last term in ( |A.9| ) vanishes, leav- 
ing only the first two. Thus, for any choice of light-like vector t', the holomorphic theta series 
P)M (t, t, t', b, c) in ( |A.lUp gives a modular completion of the Eichler integral $ p>M (r, t, b, c) (or 
conversely, the Eichler integral gives a modular completion of the holomorphic theta series). 

A. 4 An infinite sequence of mock theta series 



Finally, let us describe an infinite sequence of indefinite theta series similar to QA.7 ), obtained by 
applying the projection of the covariant derivative D (|4.9Q along a time-like vector t. This 
derivative acts on any Jacobi form of weight (f), Jj) by changing its modular weight to (f), f) + 1). In 
particular, it preserves the class ( A.1 ) of theta series and acts on them in the following way 



CP, A; r; b, c) = ^ (d {t) V, A + 1; r; b, c) , d [t) V(x) 
Using this result as well as the identity 



. t-{d x + 2tix) 
V2^ 



(x 



+ 7> 



V{x). 
(A.17) 

(A.18) 



where H m are the Hermite polynomials, one finds that the m-th derivative of QA.10 ) is given by 

(\ m/2 
k) £ 



(5j ) ) m e prf .(r,t,f ) 6 ) c) = ( 



1 



sgn((/c + by* ) — sgn((A; + b)\ 
H n {iV2^(k + 6) J>) r^ (fc+b)2 E(-c • (fc + |b)) . (A.19) 



x 



(*) 



Similarly, the m-th derivative of ( A. 11 ) is found to be 

{D^) m ^ P Ar,t,b,c) = r 2 m/2 (~ 1 ) k ' PS m (V§£|(* + &) 

fceA+/x+|p 

where we introduced the function 

S m (x) = d§ (sgn(x^) Erfc {yft\x+\) 

ml (2vr) m / 2 [sgn(x + )] m+1 H^{^\x + \) e~™+. 



E(-c- (fe + ib)) 



(A.20) 



(A.21) 
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The modular anomaly of ( A.2U| ) can be exposed by writing it as an Eichler integral similar to 
( |A.13| ) using the integral representation of S m 



/ioo / i 
E\--{k + b)[ 



(T - W. 



X 




1 (_E w—t \ 2 (fc+b)+ dw 

T)\ \ 2 w~f) [_i( w _r)]l/2 : 



m— 1 



1 ( _tt w-t \ 2 -iy^dw 



m even 
m odd 



(A.22) 



In this way, we find 



-m/2 



m! 



X) (~l) k - p ri (k+by+ q-t {k+b) - E(-c ■ (fc + ±6)) 



1 ( _E w -t \ ~2 (fc+fe)+ dw 

WS) I 2 iu-f/ [-i^-f)] 1 / 2 ) 




m even 



(A.23) 



TO— 1 



(_zl«^ 2 -y^d«, m odd 



The two theta series ( A.19| ) and ( A.2U| ) have anomalous modular transformations, but their differ- 
ence, provided the vector i! belongs to the boundary of the Kahler cone, transforms as a modular 
function of weight (0, | + m). 

The infinite set of derivatives of the theta series constructed above occurs in the analysis of 



instanton corrections to Darboux coordinates in §4.4 through their generating functions 



m=0 



m=0 



v_ 



Using 



2xy-y 2 



5> 

m=0 



ml 



(A.25) 



one finds that the first generating function is given by 



1 



fceA+jx+^P 

X 



(A.26) 



E( (ifc + + ^ ) g-2 (fe+6) E(-c • (fc + §6)) 



while the second generating function can be easily found from ( |A.23| ) 



dw 



-i(w - f )]V2 



E 



ly w — r 

4t 2 w — f 



x 



E (- 1 ) feP ((* + + ^ (fe+6) " E(-c ■ (k + ib)) . 



(A.27) 



fceA+^+^p 

Assuming that the generating parameter y transforms under SL(2, Z) with weight (0, —1), it follows 
that the difference (y) — Qp^iy) is a non-anomalous modular form, with the same weight as 
that of the original theta series ( |A.9[ ), i.e. (0, |). 
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B. Instanton corrections in the large volume limit 



B.l S-duality and mirror map 

To establish the modularity of the type IIA construction of D-instantons, we have to find instanton 
corrections to the classical mirror map ( |2.6p such that the Darboux coordinates on the twistor 
space expressed in terms of the type IIB fields transform as in ( f2~20l ). For this purpose, it will be 



useful to borrow some methods from the work |36| where a manifestly S-duality invariant twistorial 



construction of a QK manifold with two commuting isometries was provided. This is despite the 
fact that the Poincare series representation is not immediately applicable to our framework, as 
further discussed in appendix |C|. 

One of the main insights of the construction in |36|] was that the kernel entering the 



integral equations ( p. 38 ), which transforms in complicated way under S-duality, can be replaced by 
an invariant kernel 

K ^ s l(t^ + ll^l\^. = lt±l^.. (B.1) 
v ' ; f 2\t'-t l/t' + t'J V 2z'-z z> K ' 

It differs from the original one by a ^-independent term which can be absorbed into a redefinition of 
the coordinates parametrizing the HM moduli space ( A , £a, a. Essentially this redefinition provides 
a substantial part of quantum corrections to the mirror map. 

However, the kernel (|B.1|) is still not convenient for the study of the large volume limit because 
the additional i-independent term gives a divergent contribution in this limit (recall that the inte- 
grals over t! are localized at t! = ±i). Fortunately, there is a way to cure this problem. To this end, 
let us note that the invariant kernel ( P-l|) is not unique. In particular, two other natural choices 
suggest themselves, 

..dt' dt' 1 + it dz' T ^ , .dt' dt' I -it z dz' ,„ . 

K + (t,t') — = = , K~(t,t') — = = . B.2 

v ' ' V t'-tl + W z'-z K ' ' V t'-tl-it' z'z'-z v ; 

Both are S-duality invariant, differ from the original one by a t- independent term, and moreover 
are finite in the limit where t! — > =pi, respectively. The crucial difference between them and ( |B.1| ) 
is that the latter is invariant under the combined action of the antipodal map q : t i— > — 1/i and 
complex conjugation, whereas the kernels ( |B.2|) are mapped into each other 



dt' 

KHt,t>)^- 



r(i,Oy- (B.3) 



V 

This plays an important role in ensuring the reality conditions on the Darboux coordinates 

^§I=& W] = t£> ^W}=^\ (B.4) 

where Ui is the image of the patch Ui under the antipodal map. Therefore, the new kernels can be 
used to generalize the construction of modular invariant QK manifolds in |36j as follows: 



Let Z be the 2n + 1- complex dimensional contact manifold defined by the infinite covering 

Z = U + UU_UU U (u' mjn W+,n) U (u' m) „W- n ) , (B.5) 
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and transition functions 

H [ +0] = F cl (e A ]); ^[-0] = ^cl ( ^_ ]); H[ (m,n)^ = ^ ( ^ ^ gfl ) > (B 6) 

where F cl (X) = —K a b c ^-^r~ is an arbitrary cubic prepotential, A = (0, a) = 0, 1, . . . n — 1. Here 
U± are the usual patches around the poles o/P 1 , are two sefo of patches which are mapped to 
each other under SL(2, ^-transformations, andUo covers the rest. They must be chosen so thatUo 
is mapped to itself under the antipodal map and S-duality, whereas U m ,n are mapped to each other 
according to 

m'\ ( a c\ (m 



* m ,»] = ^ ( n , ) = ( 6 s ) ( J • ( B - 7 ) 



TTie holomorphic functions n are assumed to transform as 



22 



? [ G m,n] = G m,n, ^ ^cT^ + n0n - linear terTO5 + ( R8 ) 

where +reg. denotes equality up to terms which are regular in U^, ,. 

Given these conditions, the associated QK manifold Ai can be shown to carry an isometric 
action of SX(2,Z). The only difference with [KB is the presence of two sets of twistorial data, 
labeled by + and — , which are mapped to each other by the antipodal map. Physically, they can 
be interpreted as contributions of branes and anti-branes. The Darboux coordinates in the patch 
Uq then satisfy the following integral equations: 

e=± rn,n ^m,n 

S ] = S + t- l F*0O - tF«(Y) - J2 E'/ ^ K%t, ^G^, (B.9) 

c 1_ ™ ^ •/ Cfr, „ 



■5(5 + C A 4) d - (r 1 + *) (r'f^Y) + «F d (?)) - Cel (r'itf (r) - tfj'(F)) 



££'£ ^ [«*(*. o(i-«s wl (oO G 

e=± m,n °m,n 



£ 

m,n 



e 



where C^ n are contours on P 1 surrounding U^ n counter-clockwise and (r 2 , Y a , (°, Q { , a d ) are 
coordinates on M.. 

The key to prove the presence of S-duality is to find expressions of the coordinates (Y a , Q v a cl ) 
in terms of (r, b a , t a , c a , c a , ip) which ensure the proper transformations of the Darboux coordinates. 
Let us do this explicitly for the coordinates entering £ a . To this end, note that the derivative 



22 We refrain from writing the full non-linear transformation of G„„, since we are interested only in the linear 
approximation in this paper. The transformation property required for modular invariance at the non-linear level 
can be found in Ipql . 
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d;[o]G mn (t') transforms, like the transition function (IB. 81), with an overall factor of (c£°(t') + d) 1 . 
Using the property 



1 - Z 2 H> — (1 - z ) 

cr + d K ' 



(B.10) 



this factor can be converted into (c£°(£) + d) 1 by multiplying the kernel by ■ Put differently, 
the Darboux coordinate £ a would transform as in ( |2.2(J| ) if it were equal to 



if,', =Th" ~ r : -P^{f + izlf) 



1 - Z 2 



— v' 

f 2vri^ 



dz 7 1 — 2 



/2 



"-7? 



z' — z 1 — z 



— 8- C + 

2 £.a m,n 



dz' 1 — z' 2 Z 3 



( - n 



z f — z 1 — z 2 [z 



a- g: 



2 (~,'\3 m,n 



(B.H) 



It is easy to verify that this equality can be achieved if one chooses 



ra.n ~ ^m,n 



Y a = j (b a + it a ) - ' f 

m,n _ 

C3 = - (c a - nb a ) - £ 



dz 



l-z)dg ] G+ hn 



dz 



zdz 



d;[0] „ — 



, . — <^ L'-'J ^ m,n T A • Q ^£l"J ^ 771,71 

r<+ 47T1 ?<l //-<- 47nZ d 4a 



r - 07T2 

dz 



., v l - z) ^io]G m „ 



9;[0] Gv 



(B.12) 



We will leave the other relations undetermined because their derivation is more involved and they 
are not of interest for our purposes. 



B.2 Large volume limit 



We shall now use the above results, in particular ( B.2 ) and ( B.12 ), as a guidance for finding the 



mirror map between type IIA and type IIB fields. To apply them to our situation, we replace 
(C^ n ,G^ n ) by (i± 7 ,H± 7 ) where 7 = (0,p a ,q' a ,q ) with p ■ t 2 > 0. After this replacement, the 
Darboux coordinates are given by the same integral equations ( P-9|) where Q h a cl differ from 
the type IIA fields ( a , Ca, cr by the terms coming from the redefinition of the kernel. 23 Our aim here 
is to study their large volume limit t a — > 00. 

In this limit we concentrate on the infinitely small region of the P 1 fiber around the S-duality 
invariant point z = such that the combination zt a stays finite. Moreover, all integrals along l± 1 
are dominated by a saddle point at z' = for the set of charges 7 and at z' = 00 for the set of 
charges —7. Using this fact and the relations ( |B.12| ) (with the above replacements), the Darboux 
coordinates of the type IIA construction in the one-instanton, large volume approximation are given 



23 cj c1 also contains an additional classical term ^-K a t c b a b b b c 
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by Q4.2|) where the instanton contributions read as 



8$L = K'a + E 

<?A 

K = K + E 



dz' 



{z>) 



H 



A2 -7 



<1 :' ( + ^ K abc t a t b V c (z + 2z) + ~ k^ * 6 ^ ) // 



1 /" dz' 

-- jf (2i^'6 a + t a (l + 2^')) if_ 

, , ( h + q' r + q' a {rb a -c a + 2r 2 t a z') 



(B.13) 



5«' = -- Ua + c A 5a +E 



'/A 



z — z 



-\ K abc p c (ic a t b + 2r 2 & a t 6 + U(n - 3ir 2 ) + ^ (n - 3ir 2 )J # 7 

+ 7 / 7^ ((2ic a - Ar 2 b a )z' + t a (n(l + 2zz') + ir 2 (3 + 2zz'))) H_ 



Here S(' A , 5a denote instanton corrections to the classical mirror map (2.6) to be found from the 
condition of modularity. 

To study the modular properties of QB.13 ), it is more convenient to replace the last three 
expressions by the combinations introduced in (|4.5|) . Using the covariant derivative operators 
(|4.9|), they can be shown to be given by 



5C = 2mp a J p 
5£> a = 6(' a -D a J p + 1 -K abc t b 



^ C E 



</A 



dz' H~ + 



dz' 



[z>) 



A2 7 



5 + a' = t5(' - - (8a + C A 5C) + [^ ab D a D b - 2 



(B.14) 



+ )r. K abc t b p c V / d*' (6°^ + 3ir 2 ) - c a + 2r 2 t a z') H, - [ ^- ((c a - fb a )z' + r 2 t a ) F_ 
La' = f5(' -\{5a + ( A 5(' A ) ~ [^r 2 V^ + (c a - f &')£> a ] ^ + ir 2 z(p • t 2 ) E / 7^2 #-7 

<3A 

1 /" f dz' 

+ 7T7 E / dz ' ( 6V - ° a + r 2 tV ) ^7 " / TTvT ((c a + (3ir 2 - n)&V + 2r 2 t a ) H 

Zl 7. h~, J I--, \ z ) 



where we used the notation (O). All the terms independent on the P 1 variable z can be canceled 
by choosing appropriately the mirror map. This fixes the instanton corrections to the mirror map 
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in the large volume, one-instanton approximation to be 



5z a = - 
5C=0, 

K=- 



TTr. 



'/A 



dz (1 — z) dmH^ 



dz 

— (1 - z) d-io]H_ 

7. A i a 



5Co 
5a 



QA 

q A Je i 

Kabct b Re / dz (( r i - 4ir 2) &a - ( ir i + ^ 2 )t a z) d L H T 

71 QA 



(B.15) 



Plugging these relations into ( |B.14| ), one arrives at the expressions ( |4.7| ) given in the main text. 
Note that the integrals over £ 7 are Gaussian and can be easily computed. It is worth stressing that 
5z a gives a correction to the mirror map between the type IIA complex structure moduli z a and 
the type IIB Kahler moduli b a + it a . Such corrections do not arise for Dl-D(-l)-instantons, but are 



in general necessary in the presence of D3-instantons [36 



C. Relation to the Poincare series construction 



In this section, we make a first attempt at recasting the type IIA twistorial construction into the 
manifestly S-duality invariant framework developed in [[36]. The basic idea is to represent the 
partition function (|3.19|) as a Poincare series, and assign each term in the sum over cosets to a 
different transition function G m , n . 

For this purpose, recall that the partition function of MSW invariants (|3.21 ) is a vector-valued 
holomorphic modular form of negative weight f) = —62/2 — 1, therefore the positive frequency 
Fourier coefficients can be expressed in terms of the polar (negative frequency) coefficients by the 
Rademacher formula J75) 

i-t> 

1»' 



( 4:71 



0<<&<r-x(X>)/24 c=l 

where K(q , q' Q , c) is the Kloosterman sum 

-c<d<0,(c,d)=l \ / \ / 

Using this representation one can express the elliptic genus ( |3.19| ) as a Poincare series pSfl 24 
Z(T, y) = \ ]T E E ( cr + d ^ cf + d y 112 n ™tio) (-l) kp 



A-^— V^M). (CI) 



(C.2) 



xE 



1 



-k 2 



ar 



k ■ y 



cr + d ct + d 



R 



27rig 
x ~~i \c(cr + d) 



(C.3) 



24 For simplicity, we assume that the constant term of ( 3.21 ) vanishes. Moreover, the sum over c, d should be 
regulated to 1 < c < K, \d\ < K, before sending K to infinity. 



-37- 



where the factor 



R^x) 



(C.4) 



r(i-fj) 

approaches 1 exponentially fast at Re (x) —> oo and Rf,(x) ~ x 1_f) /r(2 — f)) at x — > 0. The group 
element g = ^" ^ j runs over the coset r oo \SL(2,Z), i.e. (c,d) run over pairs of coprime integers 
and a, b are chosen such that ad — be = 1. 

Applying the same strategy to the formal theta series ( |3.25| ) on the twistor space, one 
arrives at a formal Poincare series. This gives 



H { i) 



E E E E ^ 



(C.5) 



where 7 = (p, ^, fc, g ) 



G, 



c,d;7 
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(C.6) 



and the action of g G r oo \S'L(2,Z) on £°, ^ and ^' is given by ( |2.2U ). Although the series flU.Sp 
is divergent due to the sum over Dl-brane charges k, it formally transforms like a holomorphic 
modular form of weight —1. 

The representation ( |C.5|) suggests a natural Ansatz for the holomorphic transition functions 
( p.6| ) that enter the manifestly S-duality invariant twistorial construction outlined in § B.1| , namely 



G 



± 



G 



m,n;±7) 



(C.7) 



and choose the contours C^ n as the geodesic circles joining the two essential singularities of ( |C.6| ), 
namely the two roots n of the quadratic polynomial t(m^°{t)+n). In particular, it is easy to check 
that the contact potential generated by these twistorial data in the large volume limit coincides 
with fl4.14|) . Unfortunately, the transition functions ( |C.6|) do not satisfy the transformation property 
( P-8|) due to the presence of the i?-factor ( |C.4j) which transforms as a period integral. As a result, 
the SL(2, Z) invariance of the construction is not guaranteed, and indeed the Darboux coordinates 
computed from ( |C.6| ) acquire a modular anomaly. This anomaly should however be cancelable by 



a similar holomorphic contact transformation as the one described in §4.4| . We hope to put this 
construction on a more solid basis in future work. 
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